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Abstract. We establish a simple relation between curvatures of the group of 
volume-preserving diffcomorphisms and the lifespan of potential solutions to the 
inviscid Burgers equation before the appearance of shocks. We show that shock 
formation corresponds to a focal point of the group of volume-preserving dif- 
feomorphisms regarded as a submanifold of the full diffeomorphism group and, 
consequently, to a conjugate point along a geodesic in the Wasserstein space of 
densities. This establishes an intrinsic connection between ideal Euler hydrody- 
namics (via Arnold's approach), shock formation in the multidimensional Burgers 
equation and the Wasserstein geometry of the space of densities. 



In the famous 1966 paper pQ V. Arnold described the dynamics of an ideal fluid 
as a geodesic flow on the group of volume-preserving diffeomorphisms of a fixed 
domain equipped with the metric defined by the kinetic energy. He also showed how 
sectional curvature of this group enters the problem of Lagrangian stability of Euler- 
ian fluid motions. In this paper we are concerned with the exterior geometry of the 
group of volume-preserving diffeomorphisms considered as an infinite-dimensional 
submanifold of the group of all diffeomorphisms. We study its second fundamental 
form and introduce the associated shape operator. The focal points of this infinite- 
dimensional submanifold turn out to correspond to shocks forming in potential solu- 
tions of the inviscid Burgers equation on the underlying finite-dimensional domain. 
This provides a common geometric framework for both Eulerian hydrodynamics of 
ideal fluids and the phenomenon of shocks of the inviscid Burgers equation. 

More precisely, let M be a compact n-dimensional Riemannian manifold. Consider 
the group V s (M) of Sobolev class diffeomorphisms of M along with its subgroup 
T> S (M) of diffeomorphisms preserving the Riemannian volume form /i. (As usual, 
ifs>n/2 + l both groups can be considered as smooth Hilbert manifolds.) For a 
curve r](t) in T> S (M) defined on an interval [0, a] its L 2 -energy is given by 
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where the norm ||X||| 2 = (X,X)l-2 is defined by the L 2 -inner product 
(0.2) (X,Y) L 2 = [ (X(x),Y(x))dx 

J M 

on each tangent space T V T> S (M). The corresponding (weak) Riemannian metric 
on V s (M) is right-invariant when restricted to the subgroup D S JM) of volume- 
preserving diffeomorphisms, although it is not right-invariant on the whole of T> S (M). 

Arnold [1] proved that the Euler equation of an ideal incompressible fluid occu- 
pying the manifold M, 

d t u + V u u = — Vp 

for a divergence-free field u, corresponds to the equation of the geodesic flow of the 
above right-invariant metric on the group T) S JM). He also showed that Lagrangian 
instability of such flows, regarded as geodesic deviation, can be estimated in terms 
of sectional curvatures of the group of volume-preserving diffeomorphisms T> S (M) 
and provided first such curvature estimates of this group. 

Below, by regarding the group of volume-preserving diffeomorphisms as a sub- 
group of the group of all diffeomorphisms (cf. e.g. jTj), we describe its sectional 
curvatures by means of the second fundamental form of the embedding T) S (M) C 
V s (M) and relate its principal curvatures to the distance from its nearest focal 
point. Recall that the distance from a submanifold N to the first focal point along 
a geodesic in the direction normal to this submanifold gives a lower bound on the 
principal curvature radius of N (in the subspace containing the normal direction): 

dist (N, nearest focal point) = min | curvature radius of N | 

= 1/ max | principal curvature of N I . 

Hence a lower bound on the distance to a focal point provides an upper bound for 
(principal) curvatures of N. This motivates our study of the Riemannian geometry 
of the embedding T> S (M) C T> S (M) as the following sequence of implications, which 
we make precise below: 

i) Geodesies in the full diffeomorphism group T> S (M) with respect to the above 
L 2 -metric are described by solutions of the inviscid Burgers equation. The geodesies 
normal to the submanifold T) S (M) C V s (M) of volume-preserving diffeomorphisms 
are given by potential Burgers solutions, see Figure 1. 

ii) The first focal point along a (potential) Burgers solution determines the mo- 
ment when the shock wave develops. The geometry of the initial profile allows one 
to precisely estimate the lifespan of this solution. 

Hi) Focal points along normal geodesies are in one-to-one correspondence with 
conjugate points along the projection of these geodesies to the space of densities 
V(M) equipped with the Wasserstein L 2 -metric. 

iv) On the other hand, sectional curvatures of T) S (M) can be explicitly com- 
puted from the Gauss-Codazzi equations using the second fundamental form (or, 
the corresponding shape operator) of the embedding T> S (M) C T> S (M). 
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Figure 1. Diffeomorphism group T> S (M) projects to the space of 
densities V(M) with the fiber T) 8 (M)\ normals to T> S JM) are hori- 
zontal geodesies, and focal points along them correspond to conjugate 
points along geodesies in V{M). 



v) The location of the first focal point along a normal geodesic can be estimated 
in terms of the spectral radius of the shape operator of T) S AM). 

In the next three sections we describe the geometry behind relations i) — in), 
respectively. The second fundamental form, the shape operator and its spectral 
properties in items iv) — v) are described in Section HI A detailed analysis of the 
shape operator will be postponed to a future publication. In the last section we recall 
the properties of asymptotic directions on the diffeomorphism groups discussed in 
[9] within the framework of the differential geometry of T> S (M). Finally, we mention 
that the connection between the Lagrangian instability of ideal fluids and negativity 
of curvatures of the corresponding diffeomorphism group was precised by V.Arnold 
in [H, [2] , and further explored in [10] , [H] . 



1. Burgers potential solutions 

The geometric characterization of the onset of shock waves in terms of focal points 
of the embedding V s (M) C T> s (M) is based on the following observation. 
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Proposition 1. Geodesies in the group T> S (M) with respect to the L 2 -metric (EHf 
correspond to solutions of the Burgers equation: each particle moves with constant 
velocity along a geodesic in M. Geodesies normal to the submanifold T> S (M) have 
potential initial conditions. 

Proof. First assume for simplicity that M is equipped with a flat metric. Then the 
statement that particles move with constant velocity along their own geodesies in 
M follows from the "flatness" of the L 2 -metric ( 10. 2p . To see that the corresponding 
velocity field satisfies the Burgers equation denote the flow of fluid particles by 
(t, x) — > r)(t, x) and let u be its velocity field 

drj 

—(t,x)=u(t,r)(t,x)), rj(0,x)=x. 
The chain rule immediately gives 

■jpfa x ) = (d t u + Du ■ u)(t,rj(t,x)) 
and hence the Burgers equation on M 

d t u + Du ■ u = 

is equivalent to 

the equation of freely flying non-interacting particles in M. 

In the general case, any Riemannian metric on M induces a unique Levi-Civita 
L 2 -connection V on D s (M), which is determined pointwise by the Riemannian con- 
nection V on the manifold M itself (see for example [7] or [12] ). Then the same 
chain rule leads to the Burgers equation 

d t u + V u u = 

equivalent to the Riemannian version of freely flying particles: V^rj = 0. The 
latter equation has the explicit form fj + ^ij k F'j k (r])f]if] k -£j = 0, where T l - k are the 
Christoffel symbols of V in a local coordinate system x 1 , . . . , x n on M. 

Finally, observe that the tangent space to the subgroup Z>*(M) at the identity 
consists of divergence-free vector fields and hence the space of normals is given by 
gradients of H s+1 functions on M. Thus horizontal geodesies are the ones whose 
initial velocities are gradient fields u\ t= o = grad0 o . □ 



2. Burgers shock waves 

The following result describes a connection between focal points and formation of 
shocks for the Burgers equation. 
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Figure 2. The velocity profile for a solution of the one- dimensional 
Burgers equation. An inflection point of the initial profile u generates 
the shock wave. The perturbed profile u leads to a later formation of 
the shock wave. 



Theorem 2. The first focal point to the submanifold V S ^(M) C V S (M) in the di- 
rection grad O is given by the moment of the shock wave formation for the Burgers 
equation with this initial condition u\t=o = grad0 o . 

Example 3. In the 1-dimensional case the shock wave solutions of the Burgers 
equation appear from the inflection points of the initial velocity profile, see Figure 



2. Such points correspond to u" 



\t=o 



0, i.e. to 



for u 



Proof. Before the appearence of a shock wave any solution with smooth initial data 
remains smooth and the shock wave is the first moment of non-smoothness. Fur- 
thermore, for a given potential <ft , any such solution is given by a family of dif- 
feomorphisms rj(t) : x i— > exp(t grad0 o (x)) parameterized by time t. The loss of 
smoothness occurs when the potential —t <p ceases to be c-concave with respect to 
the distance function c(x,y) = d 2 (x,y)/2 on the manifold. In this case i](t) is no 
longer a diffeomorphism. On the other hand, as long as the potential is c-concave, 
the curve r](t) remains the shortest curve joining the diffeomorphisms f](0) and //(£), 
see [11], and focal points on r](t) cannot appear. 

Finally, we note that the moment when a shock wave appears corresponds to a 
focal point (i.e., to a caustic), rather than to a cut (or Maxwell) point. This follows 
from the fact that in any shock wave the first collisions start between neighboring 
particles, see e.g. [21 E]. To see why such collisions lead to the existence of a neigh- 
boring shorter geodesic, one can employ the following heuristic argument. Consider 
a shock wave generated by a solution to the 1-dimensional Burgers equation with 
the initial velocity u = grad0 o , see Figure 2. Assume that this shock wave emerged 
at time t\ and consider it at time £2 > t\- Consider a perturbation u of the original 
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solution which at the initial moment differs from the initial velocity u on a small 
segment. This segment consists of exactly those particles which are glued together 
in the shock wave between the times [ti , 1 2] ■ Perturbations u, whose velocity profile 
on the (x, w)-plane (see Figure 2) lies below that of u correspond to smooth per- 
turbations of the initial velocities of these particles, slowing them down to ensure 
that the shock wave develops after time t\ but no later than t?,. Note that such 
a perturbation defines another solution which coincides with the original one after 
time t<i but represents a shorter geodesic. Indeed, since we decreased the velocity 
vectors of the particles, the length of the new geodesic, being the integral of the 
velocity squares, can only get smaller. By arbitrarily shortening the segment 1 2 ] 
this implies that the first moment t\ of the original shock wave for the initial con- 
dition u indeed defines a focal point, rather than a cut point, since there exists a 
shorter geodesic nearby. □ 

Remark 4. In higher dimensions the shock waves are first generated from the special 
points of the initial potential u\ t =o = grad <p : singularities of type A 3 modulo certain 
linear and quadratic terms in local orthogonal coordinate charts, see [5J. The list 
of initial singularities, possible bifurcations of the shock waves and other related 
questions for the inviscid Burgers equations, can be found in [3], HI [8]. A given 
shock wave for the Burgers equation forms at a caustic corresponding to a deeper 
degeneration of this potential, see e.g. [3]. 

The most degenerate caustic points correspond to the endpoints of the spectrum 
for the shape operator Sh gra d(f> , which we describe below. The lifespan of a solution 
before the shock formation provides an estimate for the spectrum of this operator: a 
lower bound for the lifespan, and hence for the distance to the first focal point gives 
the upper bound for the (principal) curvatures of the submanifold T> S ^(M) C T> S (M). 



3. Focal points in diffeomorphisms and conjugate points in densities 

Next we would like to show that conjugate points along geodesies in the Wasser- 
stein space of densities are in 1-1 correspondence with focal points on the group of 
diffeomorphisms. Namely, consider the fibration ir : V s (M) — > V(M), where tt is the 
projection of diffeomorphisms onto the space of (normalized) densities on M. Thus, 
two diffeomorphisms T]i and r\2 belong to the same fiber if and only if 77! = i] 2 o ip 
for some volume-preserving diffeomorphism <p. (Moreover, this projection can be 
defined for more general maps by tracing how they transport the given standard 
density form on M.) The map 7r is a Riemannian submersion [1_3] and the above 
L 2 -metric on T> S (M) induces the Wasserstein (or, Kantorovich-Rubinstein) metric 
in the space of densities V(M). 

The following general observation holds. Consider a projection 71 : T> —>■ V be- 
tween two (possibly weak) Riemannian manifolds which is a Riemannian submersion 
and where T> an infinite-dimensional Lie group. Let 7 be a geodesic starting at a 
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point 7 in V and let rj be a horizontal lift of this geodesic to T>, i.e. a horizontal geo- 
desic in T> whose initial point lies in the fiber 7r _1 (7 ) and whose projection to V is 7. 
Assume, in addition, that the metric restricted to 7r _1 (7 ) is right invariant. (In our 
case, the fiber 7r _1 (7 ) is the group V s (M) of volume-preserving diffeomorphisms 
for 7 ). 

Proposition 5. The points along the geodesic 7CP conjugate to the initial position 
7 are in one-to-one correspondence with the focal points of the fiber = 7t~ 1 ( / -f ), 
regarded as a submanifold in T>, along a horizontal geodesic rj in T>. Moreover, 
the multiplicities of the conjugate points in V coincide with the multiplicities of the 
corresponding focal points in T>. 

Recall that conjugate points and their generalizations, the focal points, arise as 
singularities of the Riemannian exponential map. In a general infinite-dimensional 
space two types of such points, called epi-conjugate (epi-focal) and mono-conjugate 
(mono-focal) points, can be found. Roughly speaking, geometric significance of the 
former has to do with covering properties of the exponential map, while the latter 
are responsible for the minimizing properties of geodesies. 

Proof. The statement follows from the submersion property of ir. Since geodesies 
in the base manifold have unique lifts to horizontal geodesies in T> we can identify 
the geodesies emanating from 7 that are close to 7 and which come together near 
a conjugate point in V with the horizontal geodesies emanating from 7r _1 (7 ), that 
are close to and come together near the focal point in T>. 

Furthermore, the argument works for both mono- and epi-conjugate (or focal) 
points, since the exponential map is nondegenerate along the fibers. Indeed, shifts 
along the fibers correspond to right multiplication by an element of 7r _1 (7 ), while 
the metric in T> is 7r~ 1 (7 )-invariant. □ 

Remark 6. Geodesies in the Wasserstein space V(M) of (smooth) densities are 
projections of the horizontal geodesies in the diffeomorphism group TJ S (M). For 
a given potential — <p a they are the diffeomorphisms r](t) : x i— > exp(t grad0 o (a;)) 
parameterized by t. According to the theorem on polar decomposition on manifolds 
[TTj every non-degenerate map 77 G T> s (M) has a unique decomposition rj = gr o 
vp into a "gradient map" gr(x) := exp :r (grad0 o ) for a c-concave potential — (p Q 
and a volume-preserving map vp. Consequently, the projection of the family of 
diffeomorphisms : 1 h exp(t grad0 o (x)) remains a shortest geodesic in V (in 
the space of non-degenerate maps) as long as the potential —t 4> remains c-concave. 



4. The shape operator and focal points 

We next turn to the shape operator of the embedding T> S ^(M) C T> S (M). The 
eigenvalues of this operator are traditionally called the principal curvatures and the 
corresponding eigenvectors of unit norm are the principal directions. In our case, at 
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each point in J) S JM) there is a family of shape operators parametrized by vectors 
normal to the submanifold. 

Definition 7. The shape operator of the submanifold T> S (M) at the identity is the 
operator Sh gra _ d ^, on T id £>*(M) defined by 

(Sh gradi> (w),v) L 2 = (S id (w,v),giadip) L 2, 

where ip G H S+1 (M), w and v are in T- ld D s ^(M) and Sid is the second fundamental 
form of D S (M). By right invariance we similarly define the shape operator at any 
point rj in V S JM). 

Remark 8. Recall that the (weak) Riemannian metric on V s (M) induces on both 
the group V S (M) and its subgroup V S ^(M) the unique Levi-Civita connections V and 
V = -Pr,V, where P v := R v o P ld o R v -i is the Hodge projection P id onto divergence- 
free vector fields on M conjugated with the right translation by r] (see j7j for more 
details). 

The second fundamental form S of the submanifold T>*(M) C T> s (M) is the 
difference of the two connections, which at the identity is given by 

S id (w, v) = Q ld V w v , 

where Qi d (w) := w — Pi d (w) = grad(A _1 div w) is the Hodge projection onto the 
gradient fields on M. The second fundamental form, as well as the curvature ten- 
sors of the connections on V s (M) and D^(M), are bounded multilinear operators 
satisfying the Gauss-Codazzi equations, see [T2] , 

The above gives an explicit formula for the shape operator at the identity 

Sh gTa , dl p(w) = -P id V w gradV> . 

One can check that the shape operator is bounded and self-adjoint. 

The next result provides an estimate of the distance to the first focal point along 
a normal geodesic in terms of the spectral radius of Sh^^^. 

Theorem 9. Suppose that M is a compact manifold of non-negative sectional cur- 
vature. Let r](t) be the geodesic normal to V S ^(M) at r\ with fj(0) = grad0 o o r\ and 
let A = 1 1 5" /igrad^> 1 1 i (w) - Then rj(t) contains a focal point in the interval < t < 1/A. 

In this case our strategy to show existence of focal points along geodesies normal 
to the fibre T> S ^(M) is to obtain estimates on the second variation of the L 2 -energy 
(10. ip . The latter is given by the following formula 

E"( V )(W,W) = J* (||V^||£ a - (R(W, 77)77, W) L 2^j dt - (Sh gr!id ^ o (w ),w ) L 2 

for any vector field W along 77 such that W(0) = w or] is tangent to the fibre P*(M) 
and W(a) = 0. Here R is the curvature of the connection V on the group T> S (M). 
We discuss this theorem and its proof in detail in a future publication. 
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We point out that the conclusion of Theorem [9] holds provided that the geodesic rj 
is defined on the interval [0, 1/A]. Furthermore, this geodesic does not minimize the 
L 2 -distance between T) S {M) and any rj(t) with t > 1/A and hence the focal point is 
in fact mono-focal. 

Another corollary of the same formula for the second variation is as follows. 

Proposition 10. Suppose now that M has non-positive sectional curvature. Con- 
sider any initial vector grad0 o with a convex function <f> Q . Then there can be no focal 
points on the corresponding geodesic. 

Indeed, for such a manifold M and any W one obtains 



In particular, if <p is convex then its Hessian (V^ o grad0 o , w Q ) is non-negative and 
so E"(i])(W, W) > 0, implying the absence of focal points between r/(0) and any 
other r)(a). 



Asymptotic directions for the group of volume-preserving diffeomorhisms also ap- 
pear naturally in the context of the exterior geometry of this group. Recall, that a 
vector tangent to a Riemannian submanifold is asymptotic if the (vector- valued) sec- 
ond fundamental form evaluated on it is zero. The geodesies issued in the direction of 
this vector, one in the submanifold and the other in the ambient manifold, have a sec- 
ond order of tangency. (Note that in general two geodesies with a common tangent 
will have only a simple, i.e. first order, tangency.) Asymptotic vectors for the sub- 
manifold of volume-preserving diffeomorhisms among all diffeomorhisms are given at 
the tangent space to the identity by vector fields X satisfying div VjfA = div X = 0, 
see [5J [9] . Such fields rarely exist. For instance, for diffeomorhism groups of surfaces 
one has the following sufficient condition for non-existence of such directions. 

Theorem 11. [9j If M is a compact closed surface of nowhere zero curvature, then 
£> M (M) has no asymptotic directions. 

Proof. Consider the square length function / := g(X,X) on M. At its maximum 
point xq we get from df(x ) = that the Jacobi matrix DX is degenerate at x . 
This implies that tr(DX) 2 (xo) = — 2det [.DX(x )] = for a divergence-free field X. 
Now employing the identity 



where r stands for the Ricci curvature of the metric g, which holds for any divergence- 
free vector field X on M, we find that 




5. Asymptotic directions 



divV x A = r{X, X) + tr(M) 2 



div V x X(x ) = K(x )g{X,X)(x ) , 
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since in two dimensions the Ricci curvature r and the Gaussian curvature K coincide. 
However, for an asymptotic field X this implies that 

= K(x o )g(X,X)(x o ) 

contradicting the assumption K ^ on M, and in particular at the point Xq. □ 

Remark 12. Vanishing of the second fundamental form on asymptotic vectors im- 
plies vanishing of the projections of this form to any gradient direction. In particular, 
such vectors must be asymptotic simultaneously for all the projections of the second 
fundamental form described above. This implies the following sufficient condition 
for the non-existence of asymptotic directions: if at least in one of the gradient pro- 
jections the second fundamental form is sign-definite, then there are no asymptotic 
directions. The latter sufficient condition can also be given in terms of gradient so- 
lutions of the Burgers equation: in such a gradient direction all shocks develop only 
as t changes in the one direction (say, increases), and do not develop as t changes 
in the other direction. 
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